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Abstract 



We discuss the analytic calculation of quasi-normal modes of various types of perturbations of 
black holes both in asymptotically flat and anti-de Sitter spaces. We obtain asymptotic expressions 
and also show how corrections can be calculated perturbatively. We pay special attention to low- 
frequency modes in anti-de Sitter space because they govern the hydrodynamic properties of a 
gauge theory fluid according to the AdS/CFT correspondence. The latter may have experimental 
CN \ consequencies for the quark-gluon plasma formed in heavy ion collisions. 
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1 Introduction 

To many practitioners of quantum gravity the black hole plays the role of a soliton, a 
non-perturbative field configuration that is added to the spectrum of particle-like objects 
only after the basic equations of their theory have been put down, much like what is 
done in gauge theories of elementary particles, where Yang-Mills equations with small 
coupling constants determine the small-distance structure, and solitons and instantons 
govern the large-distance behavior. 

Such an attitude however is probably not correct in quantum gravity. The coupling 
constant increases with decreasing distance scale which implies that the smaller the 
distance scale, the stronger the influences of "solitons". At the Planck scale it may well 
be impossible to disentangle black holes from elementary particles. 

- G. 't Hooft 

Quasi-normal modes (QNMs) describe small perturbations of a black hole which is a thermody- 
namical system whose (Hawking) temperature and entropy are given in terms of its global charac- 
teristics (total mass, charge and angular momentum). They are obtained by solving a wave equation 
for small fluctuations subject to the conditions that the flux be ingoing at the horizon and outgoing 
at asymptotic infinity. These boundary conditions in general lead to a discrete spectrum of complex 
frequencies whose imaginary part determines the decay time of the small fluctuations 

3(0 = - (1) 

T 

There is a vast literature on quasi-normal modes and we make no attempt to review it. Instead, we 
concentrate on obtaining analytic expressions for quasi-normal modes of various black hole pertur- 
bations of interest. One can rarely obtain analytic expressions in closed form. Instead, we discuss 
techniques which allow one to calculate the spectrum perturbatively starting with an asymptotic 
regime (e.g., high or low overtones). In asymptotically flat space, we discuss the cases of four- 
dimensional Schwarzschild and Kerr black holes. Generalization to higher-dimensional spacetimes 
does not present substantially new calculational challenges. However, we should point out that the 
case of a rotating black hole is considerably harder than the Schwarzschild case. 

We also discuss asymptotically AdS spaces and obtain the spectrum as a perturbative expansion 
around high overtones. At leading order the frequencies are proportional to the radius of the horizon. 
When expanding around low overtones, one in general obtains an additional frequency which is 
inversely proportional to the horizon radius. Thus for large black holes there is a gap between the 
lowest frequency and the rest of the spectrum of quasi-normal modes. We pay special attention to 
the lowest frequencies because they govern the behavior of the gauge theory fluid on the boundary 
per the AdS/CFT correspondence. The latter may have experimental consequencies pertaining to 
the formation of the quark-gluon plasma in heavy ion collisions. 



2 Flat spacetime 



We start with a study of QNMs in asymptotically flat space-times. We discuss scalar perturbations 
of Schwarzschild and Kerr black holes in four dimensions. 



2.1 Schwarzschild black holes 

The metric of a Schwarzschild black hole in four dimensions is 



? , x 9 dr 2 , , , , 2GM 

ds 2 = -f{r) dt 2 + — + r 2 dQ 2 , f{r) = 1 (2) 

f{r) r 



The Hawking temperature is 

1 1 

'T' 

H ~ %kGM ~ 4nr Q 

where r$ = 2GM is the radius of the horizon. 

A spin- j perturbation of frequency 00 is governed by the radial equation 

d ( „, s d*¥\ . „, „„ 2 



(3) 



" /W d~r I /W ~dr~ ) + V (r) ^ = W W 

where V (r) is the "Regge-Wheeler" potential 

, (r)=/(r) (£<£+i> + <lzp3>) <5) 

The spin is j = 0, 1,2 for scalar, electromagnetic and gravitational perturbations, respectively. It 
is advantageous to avoid integer values of j throughout the discussion and only take the limit j — ► 
integer at the end of the calculation. 
By defining the "tortoise coordinate" 

dr ( r \ 

r + r ln --1 (6) 



f(r) ' V r o 

the wave equation may be brought into a Schrodinger-like form, 

d 2 ^> 



dr 2 



+ V(r(n))y=CQ 2x i' (7) 



to be solved along the entire real r*-axis. At both ends the potential vanishes (V — ► as r* — ► ±°°) 
therefore the solutions behave as *F ~ e ±mr *. For QNMs, we demand 

xy^ e Tia>r t ^ r# _, ±00 ( 8 ) 

assuming 9^ ft) > 0. 

2.1.1 Limit £ -»• °° 

In this case it suffices to consider the potential near its maximum. Expanding around the maximum 
of the potential (V^/w) = 0) [1], 

r m a X = lr o + 0(l/l), (9) 

we obtain 

V [r(r»)] ~Cl 2 -p 2 (r*-r*(r max )) 2 , (10) 



where 

a 2 = ^(e+ l -)r 2 o + 0(i/e) , p = ^-[< + -)+o(\/(). (id 



The solutions to the wave equation are 

x i , n=H n ( v r if3x)e'^' 2 , n = 0,1,2,. .. (12) 

where H n are Hermite polynomials. The corresponding eigenvalues are 



(O n = ^ 7 U+\ + i{n + \)\ + 0{\/t) (13) 

This result is in agreement with the standard WKB approach [2]. 

2.1.2 Limit n — ► °o 

The asymptotic form of QNMs for large n is 

^ = (2«+l)7T/ + ln3 (14) 

independent of the angular momentum quantum number £. This form was first derived numeri- 
cally [3-7] and subsequently confirmed analytically [8]. The large imaginary part of the frequency 
(Zco n ) makes the numerical analysis cumbersome but is easy to understand because the spacing of 
frequencies is 2%iTji which is the same as the spacing of poles of a thermal Green function on the 
Schwarzschild black hole background. On the other hand the real part (?R(O n ) is small. Its analytical 
value was first proposed by Hod [9]. 

The analytical derivation of the asymptotic form (14) of QNMs by Motl and Neitzke [8] of- 
fered a new surprise because it heavily relied on the black hole singularity. It is intriguing that the 
unobservable region beyond the horizon influences the behavior of physical quantities. 

We shall calculate the asymptotic formula for QNMs including first-order corrections [10] by 
solving the wave equation perturbatively for arbitrary spin of the wave. We shall obtain agreement 
with results from numerical analysis for gravitational and scalar waves [5, 11] and WKB analysis 
for gravitational waves [12]. 

Let 

y = e- ia>r *f(r*). (15) 

We have /(r*) ~ 1 as r* — > +oo and near the horizon, /(r*) ~ e 2mr * (as r* — ► — oo). Let us continue 
r analytically into the complex plane and define the boundary condition at the horizon in terms of 
the monodromy of /(r* (r)) around the singular point r = ro, 

^(r ) = e - A7t<oro (16) 

along a contour running counterclockwise. We may deform the contour in the complex r-plane so 
that it either lies beyond the horizon (Rer < ro) or at infinity (r — > oo). The monodromy only gets a 
contribution from the segment lying beyond the horizon. 
It is convenient to change variables to 

z=CO(n-inr Q ) =a>(r + r ln(l-r/r )) (17) 



(where we chose a branch such that z — ► as r — ► 0). The potential can be written as a series in y/z, 



(Q 2 ( . 2 , 3£(£+\) + l-j 2 l_2z_ 
4 Z 2 I J + 3 V (Or 



Vto = -T3l 1 -f+— 5 -\/-^- + -J ( 18 ) 



which is a formal expansion in XlyfJo. 

Now deform the contour defining the monodromy so that it gets mapped onto the real axis in 
the z-plane. Near the singularity z = 0, 

CO o 

z^- — r 2 (19) 

2r 

Choose a contour in the r-plane so that near r = 0, the positive and negative real axes in the z-plane 

are mapped onto 

argw 3k argfti 
argr = 7T — , argr= — — (20) 

in the r-plane, respectively. These segments form a n/2 angle (independent of arg ft)). 

To avoid the r = singularity, go around an arc of angle 37T/2 which corresponds to an angle 
of 3n around z, = in the z-plane. 

Considering the black hole singularity (r = 0), we note that there are two solutions, 

f ± (r)=r l± Jz ± (r) (21) 

where Z± are analytic functions of r. Going around an arc of angle of 3n/2, we obtain 

/± (^/2 r)= ^(l ±i >72 /±(r) (22) 

which is an exact result. 

To proceed further, we need to relate the behavior of the wavefunction near the black hole 
singularity to its behavior at large r in the complex r-plane. To this end, we shall solve the wave 
equation perturbatively, thus writing the wavefunction as a perturbation series in \/\fco. 

At zeroth order, the wave equation reads 



^2^(0) 

~lz>~ 
Two linearly independent solutions are 



+ 0-jir + l ) * (0) = ° (23) 



fi°\z)=e^f=e^^ j/2 (z) (24) 

in terms of Bessel functions. We deduce the behavior at infinity (z — > °°) 

/i 0) (z)~e fe cos( z -7T(l±j)/4) (25) 

The boundary conditions imply f(z) ~ const, as z —> °° along the positive real axis in the z-plane. 
Therefore, we ought to adopt the linear combination 

f®=jP-e-*Wj®>~*VzHU{z) (26) 

(in terms of a Hankel function). As z —> °°, we obtain 

f(0) ( z ) „ _ e -n(\+j)ilA sinfaj/2) (27) 



a constant, as desired. 

Going along the 3 % arc around z = in the z-plane, we have 

f(0)/,3»,l - -3x(l+j)i/2 f f(0) M _ _-7wji/2 f (0) 



fW (e in, z) = e^ + M z (/}"' (z) - <T ' W / w (z) J (28) 

As z — ► oo, 

/(0) ( z ) „ e -*Q+W sm(3nj/2) + g^ 1 --/)'/ 4 si n(2^> 2 ' z (29) 

The monodromy to zeroth order is 

•*(*>) = - Sin{ f nj/ ^ = -(1 +2cos(;rj)) (30) 

sm(Kj/2) 

leading to a discrete set of complex frequencies (QNMs) [8] 

— = (2n+ l)jri + ln(l +2cos(nj)) + ffiX/y/n) (31) 

Next, we calculate the first-order correction to the above expression [10]. Expanding the wavef unc- 
tion in l/y/a, 

<p = y(°)+ l *y (1) + ^(l/co) (32) 

V-C0ro 

the first-order correction obeys 



^ + fc£ + .V>=^«V*«> 



(33) 



where 

5V(z) = ^ + ^V[r(z)] (34) 

Two linearly independent solutions are 

yW (z) = <^f (z) f ^i 0) 5y*i 0) - #>p(_ 0) (z) f ^f SV^f (35) 



where *€ = l?, and the integral is along the positive real axis on the z-plane (z > 0). We obtain 



n(jij/2) 

the large-z behavior 



/(i)/ 



^ (z) ~ c_ ± cos(z - tt(1 + j)/4) - c +± cos(z -n{\- j)/A) (36) 

where 



c±± = <tf [ vfSVWf (37) 

Jo 



To obtain the small-z behavior, expand 

s ^.. ^-^ ,^.1 (38) 

6v-2fi)ro 

It follows that 

*i 1) =z 1 ^' /2 G ± (z) + ^(l/a)) (39) 

where G± are even analytic functions of z. 



For the desired behavior as z — ► °°, define 

T - Y" ' -r ' ■ { T: ' ' - e -^/2»J/(l) + g -*ji/2^(0) | + (4Q) 



where <^ ~ ^(1) and dots represent terms of order higher than 0{\/^/~go). By demanding *P ~ e !z 

as z — > +00, we fix 

§ = £+ + £_ , £+ = c ++ e*# 2 - c+ _ , £_ = c __e-^/2 _ c+ _ (4i) 

Then the requirement f{z) = e iZX P(z) ~ const, as z — ► °° yields 

-7r(l+y)«74 r :_/'_../oN J 1 €- 



/(z) ~ -e-*V + »«* sm(xj/2) 1 - -^= (42) 

In the neighborhood of the black hole singularity (around z = 0), going around a 3n arc, we obtain 

^\e 3ni z) = t?*V±M 2 yV> ( Z ) (43) 

therefore 

¥( e 3 ™ z ) = *(°)(e 3TO z ) 

As z — > 00 along the real axis, 

f(z) ~ e-*( 1+ ^'/ 4 sin(3wy/2)(l ' 



+ ^(l-7>'/4 sin ( 27r j) (l _ 1 gL 2 fe 

where 

A = i_I e ^/2 ( ^ + + ^_ _ ^ cot(3 ^/2)) (45) 

and B is not needed for our purposes. The monodromy to this order reads 

■*(*>) = ~ ^§ { 1 + i^= e^' 2 «_ - ^ + § cot(3*y/2)) } (46) 

sm(7rj/2) [ 2^-cor J 

leading to the QNM frequencies [10] 

co e xj>/ 2 

-?■ = (2n + l);r/ + ln(l +2cos(7r;)) + . — (£. - g + + gcot(3ft//2)) 

•'ff y/n+1/2 

+ @{l/n) (47) 

which includes the @(\/y/n) correction to the ^(1) asymptotic expression (31). 
For an explicit expression, use 

,00 / / 7r72r( v+M+1/2 ) 

V ^ s / ^- 1/2 Mz)Mz) = r( _ v+ ^ + V^ v ^ 2) -^^- (48) 
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We obtain 

c±± = * 3 »fi 1 Kl~™ s(±jfr±m (49) 

12y2 sm(7rj/2) 
therefore 

e e e /„ -/«n /, ., 3^+ 1) + 1 - /' 2 sin(27T./) 

x r 2 ( 1 /4) T( 1 /4 + j/2) T( 1 /4 - ;/2) (50) 

where we also used the identity r(y)r(l — v) = t f ,■, ■ This expression has a well-defined finite 
limit as j — » integer. 

For scalar waves, let 7 — > + . We obtain 

^ = (2. + l)^ + ln3 + ^i4= ^ + ^ + 1/3 r 4 (l/4) + ^l^) (51) 

which is in agreement with numerical results [11]. 
For gravitational waves, we let j — ► 2 and obtain 

| - (2„ + lM + ta 3 + JlL. 1+JjJ r-(l/4) + ^(l/„) (52, 

which is in agreement with the results from a WKB analysis [12] as well as numerical analysis [5]. 

2.2 Kerr black holes 

Extending the above discussion to rotating (Kerr) black holes is not straightforward. Bohr's corre- 
spondence principle 

8M = hSico (53) 

and the first law of black hole mechanics 

8M = T H 8S BH + &8J (54) 

imply the asymptotic expression[9] 

3i(D = T H ln3 + mQ (55) 

where m is the azimuthal eigenvalue of the wave and Q. is the angular velocity of horizon. In 
deriving the above, we identified 8Sbh = In 3 [13]. Even though the above result has the correct limit 
as Q. —> (in agreement with the Schwarzschild expression (14)), it is in conflict with numerical 
results [14] indicating SRft) w mQ.. 

To resolve the above contradiction, we shall obtain an analytic solution to the wave (Teukolsky 
[15]) equation which will be valid for asymptotic modes bounded from above by I /a, where 

a = lA (56) 

with / being the angular momentum and M the mass of the Kerr black hole. The calculation will be 
valid for a <C 1 which includes the Schwarzschild case (a = 0) [16]. Our results will confirm Hod's 
expression (55) and not necessarily contradict numerical results (the latter may still be valid in the 

8 



asymptotic regime \/a < (0). In the Schwarzschild limit (a — ► 0) the range of frequencies extends 
to infinity and our expression reduces to the expected form (14). 
The metric of a Kerr black hole is 

, 2 / 2Mr\ j2 4Marsm 2 6 J ^ £ J 2 
ds = - I 1 — J dt H dtdtp + — dr L 



2Ma 2 rsm 2 d 



1.2 



+ Zdd z + sin z d I r z + a z + Uf (57) 

where £ = r 2 + a 2 cos 2 6 , A = r 2 — 2Mr + a 2 = (r — r_)(r — r + ) and we have set Newton's constant 
G = \. The angular velocity of the horizon and Hawking temperature, respectively, are 

fl = — , T H = 1 - r ~ /r+ (58) 

2Mr+ ' 87TM 

2.2.1 Massless perturbations 

Massless perturbations are governed by the Teukolsky wave equation [15] 

'r 2 + a 2 ) 2 , , \ d 2x i> AMar d 2x ¥ fa 2 1 \ d 2x ¥ 
—!- - a 2 sin 2 -=-=- + — — -^—r + i-r- 



dt 2 A dtdtp \A sin 2 6 J dtp 2 

1 d [ AS+l dV\ „ (Mir 2 -a 2 ) . \ d*¥ 

A s dr V dr J \ A J dt 

I d ( . n d^\ „ (a(r-M) icosd\ d*¥ , , in ,„, 
'sin0 — -2s (-^-j — - + — -5— 1 — + (s 2 cot 2 d-s)*¥ 



sin0 dd\ dd J V A sin 2 0/ 50 



(59) 



where s = 0,-1,-2 for scalar, electromagnetic and gravitational perturbations, respectively, a 
Writing the wavefunction in the form 

y = e -ia>t e »«$S(d)f(r) (60) 

we obtain the angular equation 

1 . ,., / 9 9 9 m 2ms cos 6 7 9 \ 

sine S) + a 2 co 2 cos 1 6 = 2acosco&6 = s 2 cot 2 6 )S 

sin0 v V sin 2 sin 2 / 

= -{A + s)S (61) 

where A is the separation constant (eigenvalue) and the radial equation 

l(A' +1 /)' + V(r)/= (A + a 2 (» 2 )f (62) 

where the potential is given by 

. , (r 2 + a 2 ) 2 co 2 - AaMrcom + a 2 m 2 + 2ia(r - M)ms - 2iM(r 2 - a 2 )(Os 

V(r) = \-2ir(Os (63) 

A 

Let us simplify the notation by placing the horizon at r = 1, i.e., by setting 

2M=l+a 2 , r_=a 2 , r+ = l (64) 



and solve the two wave equations by expanding in a. We shall keep terms up to 0(a) assuming ft) 
is large but bounded from above by \/a, (1 < ft) < \/a). Thus ft) is in an intermediate range which 
becomes asymptotic in the Schwarzschild limit a — ► 0. 

The solutions to the angular equation to lowest order are spin-weighted spherical harmonics 
with eigenvalue 

A = £(£+l)-s(s + l) + 0(a<o) (65) 

Near the horizon (r — ► 1), 

/(r)~(r-l) A , X = /(ft)- am) + 0(1 /ft)) (66) 

At infinity (r — ► <»), /(r) ~ e' mr . Introducing the "tortoise coordinate" 

z = (Or+(co-am)\n(r- 1) (67) 

the boundary conditions read 

f(z)~e ±iz , z ^±oo (68) 

From the boundary condition at the horizon we deduce the monodromy for the function &(z) = 
e' z f(z) (notice that & ~ const, as z, — ► +°°) around the singular point r = 1, 

jj[(\)=e An ^- a ^ + ff(a 1 ) (69) 

To express the radial equation in terms of the tortoise coordinate, define 

/M-y* ,, R(r) ■ TO 

Ao = r(r - 1) (note A = Aq + 0(a 2 )). Inverting z = z{r), 



— + 0(1/(0) (71) 

ft) 

the radial equation to lowest order in l/^fco in terms of R reads 

d 2 R ( 3 is 4-s 2 -4iams) 

to be solved along the entire real axis. This is Whittaker's equation. The solutions may be written 
as 

M K)±ll (x) = e~ x/2 x ±fl+l/2 M(- + jx - K,\ +2jx,x) (73) 

where K = ^, \i 2 = £i£+J£^Z 5 M K< ±n is Kummer's function (also called <1>) and we set x = 2iz. We 
need to introduce Whittaker's function 

W K ,» (x) = r ^ ( " 2M) - M K ,, (x) + ^ M K ,_„ (*) (74) 

r(i-/i-?c) r(±+ju-K) 

due to its clean asymptotic behavior, 

W K ^{x) ~ e~ x/2 x K (l + 0( 1/jc)) , |*H°°- ( 75 ) 
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We may compute the monodromy by deforming the contour as before. Going around an arc of 
angle 3 71, we have 

M K , ± ^(e M x) = -ie ±3 ^M_ K ^(x) (76) 

where we used M(a,b, —x) = e~ x M{b — a,b,x), therefore 

W K ,,(e™x) = -UP* _T ( ~ 2/X) , M_ rjf (*) - UT**» fM M {x) (77) 

To find the asymptotic behavior, we need 



As be I —>■ °°, we obtain 



(78) 



W K ^ (e 37ti x) ~ Ae*/ V + Be~ xl \- K (79) 



where 



,3«i/i r (-2/*) r(i+2/i) * 



.^ v r-, v . ^ g -^ + ( ^ } (g0) 



and B is not needed for our purposes. After some algebra, we deduce 

A = -(l+2cos7Ts) + ^(a 2 ) (81) 

where we used the identities T(l — x)Y(x) = ^^, r(j +x)r(^ — x) = — ^— . The monodromy 
around r = 1 is 

fl\ =e Mw-ma) =A (82) 



therefore [16] 

1 1 

9to = — - ln(l + 2cos to) + ma + <^(a 2 ) (83) 

4n 

in agreement with Hod's formula for gravitational waves (s = —2) in the small-a limit (in which 

Q. ss a, Th ~ 4^). However, it should be emphasized that these are not asymptotic values of QNMs 

but bounded from above by I /a. 

2.2.2 Massive perturbations 

The case of massive perturbations is interesting because it reveals instabilities. As is well-known, 
the Schwarzschild spacetime is stable against all kinds of perturbations, massive or massless which 
makes the Schwarzschild geometry appropriate to study astrophysical objects. On the other hand, 
Kerr spacetime is stable against massless perturbations but not against massive bosonic fields [17]. 
The instability timescale is much larger than the age of the Universe so the problem is not expected 
to have observable consequencies. Nevertheless, the study of instabilities is an important subject 
and QNMs provide an indispensable tool. 

For a massive scalar of mass /x, the radial wave equation reads 

A— + <^ — i 1 juV-aW-^+l) }R = 0, (84) 



dr\dr) 
We are interested in solving this equation for a small mass and low frequencies (;U, ft) <C 1/M) [17]. 

11 



Away from the horizon (r 3> M), we may approximate by 



<i 2 , , 



2 2MjU 2 ^ + 1) 



rR = , k 2 = n 2 -(0 2 



The solution to this equation is given in terms of a confluent hypergeometric function, 

R(r) = (2kr) l e- kr U{l+ 1 -M^l 2 /k,2{i + l),2kr) 
Near the horizon (r <C £/|fc|), we may approximate by 



Z (z+l)| 






+ [ J P 2 -^+l)z(z+l)]/? = 



(85) 



(86) 



(87) 



where P= ™-™r + '» * - I=r± 



-,z : 



function, 



- 1 -. The solution to this equation is given in terms of a hypergeometric 



*(*) 



z+1 



(P 



! + l;l-2iP;z+l) 



(88) 



Matching the two expressions in the overlap region (M <C r <C £/|fc|), we obtain the frequencies 

(Dn^ll + ijn , «GN (89) 

where 



% = «Wi(MM) 4(m: 



2 2(2i-+l) (2j , + 1+M) , (j ,, ) 2 



am 



M-2Air + j 



n 



-, / a 2 \ ( am \ 2 



(90) 



and <^„ - ^ +1+n )2(m)( 2 £ + i)2 n! ((2i)!)4- 

For m > 0, we have y„ > yielding an instability. For the fastest growing mode (with i = 1, 

m = 1 , « = 2 (2p state)) we have 

_ 1 24 

T ~7~aAi 2 (/uM) 7 ( } 

which is generally large. 

Notice that there is no instability in the Schwarzschild limit (a — ► 0) and for massless perturba- 
tions (;U — ► 0); in both cases, y ^ and therefore the lifetime x — ► °°. 



2.3 Half-integer spin 

In the case of a perturbation of half-integer spin we need to solve the Teukolsky equation [18] with 
potential 

vw=/w (^±l) + I) + ^_^ + il (92) 

\ y y I y y tX y 

where j is the spin of the perturbing field (e.g., j = 1/2 for Dirac fermion). We shall set ro = 1, for 
simplicity, so f(f) = 1 — -. 

Expanding around the black hole singularity z = (Or* = 0, 



1 



^nz) = -± 



30' 4-/ 



a> 



+ ■ 



,.c/ 



2z 16z 2 C0V2 z 3/2 



■ + 0(\/<o) , srf 



!+!) + 



1-7 2 



2^2 



(93) 
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we obtain the zeroth-order wave equation 

d 2x ¥ 



dz 2 



+ 



3i7 4-f 
2z 16z 2 



¥ = (94) 



whose solutions are the Whittaker functions 



¥°\z)=M x ^(-2iz) , 1 = ^-1, ] x = ] - (95) 

The calculation of the monodromy as before leads to the modes [18] 

^L = -(2 w +l)7n + rn(l+2cos7rj) + ^(i/Vn) (96) 

in agreement with the result for integer spin ( which came from the Regge- Wheeler equation). For 
a Dirac fermion, j = 1/2, so asymptotically, the real part vanishes. 

The first-order correction may also be calculated as before [19]. The result is 

— = -(2«+l)7n' + rn(l+2cos7T/) 
Th 

2i . b+A-B-+b-A+B + „, , s 



yZ]^ ^ g- 4 ^A + B_ - e 4 ^A_B + 

where 

£/ r°° dz 

h± = 4u.lo Jj2 Mx ^~ 2iz)Mx ^~ 2iz) (98) 

and A± = „.[, ^ e'^i^"^), S ± = — |— ^ e~ mX . This result appears to be a complicated 

1 (2±/l-t-/l) 1 \2±^l— A.) 

function of j, so let us look at specific cases. 
For 7 = 1/2 (Dirac fermions), we obtain 



T H - {2n+l)ni+ 2-Tn{ e+ 2) ^{i +W») <"> 



which is in good agreement with numerical data [19]. 
For j = 3/2, we find 



T H 
so there are no first-order corrections to the spectrum. 

For j = 5/2, we have 



%± = -(2n + l)ni + 0(l/n) (100) 



^ = -(2n+l)m+ l -±L^r'[ 1 ) + 0(\/n) (101) 



etc. 

All of the above spectra agree with the general expression we obtained for integer spin using 
the Regge-Wheeler equation. The relation of the latter to the Teukolsky equation is worth exploring 
further. 
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3 Anti-de Sitter spacetime 

According to the AdS/CFT correspondence, QNMs of AdS black holes are expected to correspond 
to perturbations of the dual Conformal Field Theory (CFT) on the boundary. The establishment 
of such a correspondence is hindered by difficulties in solving the wave equation governing the 
various types of perturbation. In three dimensions one obtains a hypergeometric equation which 
leads to explicit analytic expressions for the QNMs [20, 21]. In five dimensions one obtains a Heun 
equation and a derivation of analytic expressions for QNMs is no longer possible. On the other 
hand, numerical results exist in four, five and seven dimensions [22-24]. 

3.1 Scalar perturbations 

To find the asymptotic form of QNMs, we need to find an approximation to the wave equation valid 
in the high frequency regime. In three dimensions the resulting wave equation will be an exact 
equation (hypergeometric equation). In five dimensions, we shall turn the Heun equation into a 
hypergeometric equation which will lead to an analytic expression for the asymptotic form of QNM 
frequencies in agreement with numerical results. 

3.1.1 AdS 3 

In three dimensions the wave equation for a massless scalar field is 

w~r 3r (? ( ! - it) **) - ^5^ + M 2 * = ° (102) 

Writing the wavefunction in the form 

r 2 
<p = e i{m - px) y(y), y = 4 (103) 

the wave function becomes 

y 2 (y-l){(y-l) x ¥')' + cb 2 yV + p 2 y(y-l)V = (104) 

to be solved in the interval < y < 1, where 

coR 2 co . pR p 

co = = , p = — = . (105) 

2r 4nT H 2r 4nRT H 

For QNMs, we are interested in the solution 

V(y)=y{i-y) iA 2F l (i + i(&+p),l + i(&-py,2;y) (106) 

which vanishes at the boundary (y — ► 0). Near the horizon (y — ► 1), we obtain a mixture of ingoing 
and outgoing waves, 



*¥~A + (\-y)- im +A4\-y) 



+i& , _ n±2icb) 



- r(i±i(&+p))r(i±i(&-p)) 

Setting A_ = 0, we deduce the quasi-normal frequencies 

cb = ±p-in , « = 1,2,... (107) 

which form a discrete spectrum of complex frequencies with 3 CO < 0. 
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3.1.2 AdS 5 

Restricting attention to the case of a large black hole, the massless scalar wave equation reads 

^f(r)d^)--^ ( - ) d^-^V^ = 0, f(r) = l- r l (108) 

Writing the solution in the form 

r 2 

<£ = e Km-p-x)^{y) j y= (109) 

1) 

the radial wave equation becomes 

(v 2 -l)(v(/-l)^)'+(^/-^(3' 2 -l))^ = (HO) 

For QNMs, we are interested in the analytic solution which vanishes at the boundary and behaves as 
an ingoing wave at the horizon. The wave equation contains an additional (unphysical) singularity 
at v = — 1, at which the wavef unction behaves as *P ~ (v + 1) ±05 / 4 . Isolating the behavior of the 
wavef unction near the singularities y = ±1, 

V{y) = {y-l)- m l A {y + l) ±& l A F ± (y) (111) 

we shall obtain two sets of modes with the same 3d), but opposite 9?d). 
F± (y) satisfies the Heun equation 



v(/-l)F±" + {(3-i|id>)v 2 -^%-l}^ 



± 

+{f(4* i -'>-^ i 4-?h - ° ^ 

to be solved in a region in the complex v-plane containing \y\ > 1 which includes the physical 
regime r> 77, . 

For large d), the constant terms in the polynomial coefficients of F' and F are small compared 
with the other terms, therefore they may be dropped. The wave equation may then be approximated 
by a hypergeometric equation 

in the asymptotic limit of large frequencies a). The acceptable solution is 

F (x) = 2 Fi(a+,a-;c;(y+l)/2) , a ± = \-'^cb±\ , c = \±\cb (114) 

For proper behavior at the boundary (y — v 00), we demand that F be a polynomial, which leads to 
the condition 

a + = — n , n=\,2,... (115) 

Indeed, it implies that F is a polynomial of order n, so as y —> 00, F ~ y" ~ y~ a+ and *P ~ 
y -iw/4 y ±a/4 y -a + „ y -2 f as ex p ecte d. 
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We deduce the quasi-normal frequencies [25] 

<b = -^- = 2n(±l-i) (116) 

4nT H 

in agreement with numerical results. 

It is perhaps worth mentioning that these frequencies may also be deduced by a simple mon- 
odromy argument [25]. Considering the monodromies around the singularities, if the wavef unction 
has no singularities other than y = ±1, the contour around y = +1 may be unobstructedly deformed 
into the contour around y = — 1 , which yields 

Ji[(\)J?{-\) = \ (117) 

Since the respective monodromies are ^#(1) = e n6> ' 2 and ^#(—1) = e^ in ®' 2 , using 3ft) < 0, we 
deduce ft) = 2n(±l — i), in agreement with our result above. 

3.2 Gravitational perturbations 

Next we consider gravitational perturbations of AdS Schwarzschild black holes of arbitrary size 
in d dimensions. We shall derive analytic expressions for the asymptotic spectrum [26] including 
first-order corrections [27]. Our results will be in good agreement with numerical results. 
The metric is 

/7»- 2 »-2 Tn 

(118) 





ds 2 = -f(r)df 


l . dr 2 


-1- r 2 d0 2 , „ 


f(r) = 


r 2 

R 2 + 1 - 


2\l 




f(r) 


-\- r u\L d _2 , 


vd—3 


The radial 


wave equation can 


be cast into a Schrodinger-like form, 








d 2x ¥ 
dr 2 


+ V[r(n)]y -- 


= (0 2x ¥ , 






in terms of the tortoise coordinate defined 


by 














dr.* 1 









dr f(r) 



(119) 



(120) 



The potential V for the various types of perturbation has been found by Ishibashi and Kodama [28]. 
For tensor, vector and scalar perturbations, we obtain, respectively, 

Mr) =f(r) { eii + "- 3) + <*-y ('> + ^-^^ } (121) 

V(r\ f(r\\ ^ l + d - 3) I ( d - 2 )( d - 4 )f( r ) rf"{r) \ 

Vv(r)=/(r) \ ? + 4? 2{d-3)j (122) 
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Vs(r) 



4r 2 



£{£ + d-3)-{d-2) + 



(d-l)(d-2)n 



-2 



d{d-l) 2 {d-2yii 6{d-l){d-2) 2 {d-4)[£{£ + d-3)-{d-2)}^i 



R2 r 2d- 



g2 r d-5 



(d-4)(d-6)[£(£ + d-3)-(d-2)] 2 r 2 2{d - l) 2 (d-2)V 



R 2 



r 3d-9 



+ 



+ 



4(d-\)(d-2)(2d 2 -Ud+\S)[£(£ + d-3)-(d-2)}ix 2 

p2d—6 

{d - \) 2 {d - 2f(d -4){d- 6)jU 2 6(d -2){d- €)[£{£ + d-3)-{d- 2)] 2 jU 



r 2d— 6 



yd— 3 



6(d-\)(d-2) 2 (d-4)[£(£ + d-3)-(d-2)\n 

yd— It 



+4[£(£ + d-3)-{d-2)f + d{d- 2) [£{£ + d-3)-(d-2)f 



Near the black hole singularity (r ~ 0), 



1 



£?T 



Vt = -T2+- 

4r * [-2(d-2)n] — 



(d-3) 2 £(£ + d-3) 

— i* +..-, J&j = - — , — + ■ 



rf-1 

rf-2 



2(2^-5) 



d-2 



sfv 



4r * [-2(J-2)Ai]3=2 



and 



Vs 



d-\ 

d-2 



+ 



+ ... 



. d 2 -Sd+\3 £{£ + d-3) 
A/ = -^rm — tt^ — I - ■ 



2(2^-15) 



d-2 



d-\ 

d-2 



4r * [-2(d-2)n]** * 



+ ..., 



where 



. (2d 3 -24d 2 + 94d-U6) (d 2 -ld+l4)[£(£ + d-3)-(d-2)] 
£*S = tttt-, — ^t~, — — 1- ■ 



(123) 
(124) 

(125) 

(126) 



4(2d-5)(d-2) (d-\)(d-2) 2 

We have included only the terms which contribute to the order we are interested in. We may 
summarize the behavior of the potential near the origin by 



4r 2 



(127) 



where j = (2) for scalar and tensor (vector) perturbations. 
On the other hand, near the boundary (large r), 



V 



jl-l 



4(r*-r*) 2 



+ ... , r* 



dr 
W) 



(128) 



where _/„ = d — 1 , d — 3 and d — 5 for tensor, vector and scalar perturbations, respectively. 
After rescaling the tortoise coordinate (z = (Or*), the wave equation to first order becomes 



J^ + C0-^5^W / = 0, 



where 



JCC\ 



dz 2 



/ 2 -l 

4z 2 



j£f =-s/ z ~* J1 . 



(129) 



(130) 
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By treating J^J as a perturbation, we may expand the wave function 

¥(z)= x i'o(z) + G)-& x i'i(z) + ... (131) 

and solve the wave equation perturbatively. 
The zeroth-order wave equation, 

■Wo(z)=0, (132) 

may be solved in terms of Bessel functions, 

x i , Q(z)=A l ^J J (z)+A 2 ^.N J {z). (133) 

2 2 

For large z, it behaves as 

¥o(z) ~ J-[Aicos(z-a+)+A 2 sin(z-a + )] 

1 (Ai - iA 2 )e- ia+ e iz + -L (Aj + iA 2 )e +ia+ e~ iz 



2% V2n 

where a± = f (1± j)- 

At the boundary (r — > oo), the wavefunction ought to vanish, therefore the acceptable solution is 

¥ (r*) = B^coin - r„) Jj± (co(n - h)) (134) 

2 

Indeed, *F — > as r* —> f „, , as desired. 
Asymptotically (large z), it behaves as 

¥(r,)~\/-*cos[fi>(iv-r,)+/3], j3 = ^(l+j 00 ) (135) 

V 7T 4 

We ought to match this to the asymptotic form of the wavefunction in the vicinity of the black-hole 
singularity along the Stokes line 3z = 3 (cor*) = 0. This leads to a constraint on the coefficients 
Ai, A 2 , 

Aitan(a>r*-/3-a + )-A2 = 0. (136) 

By imposing the boundary condition at the horizon 

^{z)~e iz , z^ — , (137) 

we obtain a second constraint. To find it, we need to analytically continue the wavefunction near 
the black hole singularity (z = 0) to negative values of z. A rotation of z by — % corresponds to 
a rotation by — -£^ near tne origin in the complex r-plane. Using the known behavior of Bessel 
functions 

J v {e~ in z) = e- inv J v {z) , N v {e~ in z) = e iKV N v (z) - 2/cos 7TV/ V (z) (138) 

for z < the wavefunction changes to 

^o(z) = e-^' +1 )/ 2 ^/=z" { [A! - j(l + e inj )A 2 ] Ji{-z) +A 2 ^N i (-z)} (139) 

whose asymptotic behavior is given by 

-in(j+\)/2 . . . -in(j+\)l2 

*¥ ~ =— [Ai - j(l +2e jm )A 2 ] e~ lz + =— [A x - iA 2 ] e lz (140) 

\/2n \/2n 
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Therefore we obtain a second constraint 

A 1 -i{l+2e j7ti )A 2 = . 

The two constraints are compatible provided 

1 -i(l+2e^) 

tan(&)f* -j3 — oc+) -1 

which yields the quasi-normal frequencies [26] 



fl)r* = -(2 + y + y»)-tan ' 1+ ^,. +^ 



(141) 



(142) 



(143) 



The first-order correction to the above asymptotic expression may be found by standard pertur- 
bation theory [27]. To first order, the wave equation becomes 



Wi+,Wo=0 



The solution is 



W 



rz V?J J (z')<m X i'o(z') rz sR Nj (z')^*o(z') 

*i (z) = ^ (z) y o *' — 2__ ^ (z ) y o ^ — - 

where W = 2/n is the Wronskian. 

The wavefunction to first order reads 



where 



a\{z) 

ai{z) 

b(z) 



^G)-fc§ f Z dz'z'-^Jj(z')J;(z') 

2 Jo 2 2 

2 J() 2 2 



Z 1 



%£^ d-3 , , , 

— fiT3=! / dzV d - 2 Jj(z')Nj(z') 

2 Jo 2 2 



and ^e/ depends on the type of perturbation. 
Asymptotically, it behaves as 



¥(z) ~ J-[AiCos(z-a + )+A 2 sin(z-a + )] , 



where 



Aj = [l-fc]Ai-a 2 A 2 , A' 2 = [l+^]A 2 + aiAi 
and we introduced the notation 

d\ = ai(°°) , a 2 = a 2 (°°) , b = b(oo). 

The first constraint is modified to 

Aj tan(wr* — j3 - a+) - A 2 = 



(144) 



(145) 



¥(z) = {A! [1 - 6(z)] - A 2 a 2 (z)} ^ (z) + {A 2 [l + *(«)] +A x ai (z)} y/zN L (z) (146) 



(147) 
(148) 
(149) 

(150) 
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Explicitly, 



[(1 -£)tan(cof* -/3 — CC+) — a\]A\ - [1 +& + a 2 tan(ft)r„ ! -j3 - a + )]A 2 = (151) 



To find the second constraint to first order, we need to approach the horizon. This entails a rotation 
by —71 in the z-plane. Using 



a x (e- l% z) 



a 2 (e- in z) 



^e- inj ai{z) 



e d ~ 2 



2^J 



e ,7IJ a 2 (z) - 4cos 2 -fa, (z) - 2i(l + e in] )b{z) 



b{e- in z) = e- , '^[Kz)-*(l+«" toy )aiW] 
in the limit z — ► — °° we obtain 

¥(z) ~ -ie- ijn/2 B Y cos(-z - a+) - «^ /2 B 2 sin(-z - o+) 



(152) 



where 



B l = A x -A x e- in — 2 [b-i{\+e- iKi )a x \ 



:~r d— 3 



e + ''^a 2 - 4cos 2 ^ai - 2/(1 + e +I ' 7t - / > 



i(l +e i7Cj ) A 2 +A 2 e- in —2[b-i(\+e- i7lj )a l ] +A 1 e- i7l — e - illj a l 



B 2 



A 2 +A 2 e- i7l —^[b-i{\+e- i ^)a l )+A l e- i7l ^e- 



• Cl\ 



Therefore the second constraint to first order reads 

[1 - e~ i7C< & (ioi + b)]Ai - [i(l + 2e'^') + <T /7C ^ ((1 + g^')ai + e i7lj a 2 - ib)]A 2 = 
Compatibility of the two first-order constraints yields 

1 + b + a 2 tan(fi)f* - j3 - a+) i(l + 2<? I ' 7C ->) + «-''*^ ((1 + e^')ai + e^'a 2 - i&) 
1 - b) tan(wf* - /$ - a+) - ai 1 - e - ' 71 ^ (j^ + £) 

leading to the first-order expression for quasi-normal frequencies, 

K 



(153) 



(154) 



con 



-- {bib - 2ie~ i7l ^b- 9a x + e~ in ^a\ +a 2 -e 



l a 2 



where 



Cl\ 

a 2 

b 






~ d - 2 r (^l) r (2 + 2(d-2)) 
r ( 2(d-2) ) r (l+ 2{d-2). 



„ „ n(d-3) n( d-3 

l+2cot^h — zr-cot- -] + 



2{d-2) 2 



d-2 



Cl\ 



n(d-3) 

COt 2(d^) dl 



Thus the first-order correction is ~ <^(ra~~). 

The above analytic results are in good agreement with numerical results [29] (see ref. [27] for a 
detailed comparison). 
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3.3 Electromagnetic perturbations 

The electromagnetic potential in four dimensions is 



V EM = e -^^f(r). (155) 



Near the origin, 



; 2 -l eu+i)r* 3/2 

where 7 = 1. Therefore we have a vanishing potential to zeroth order. To calculate the QNM 
spectrum we need to include first-order corrections from the outset. Working as with gravitational 
perturbations, we obtain the QNMs 

i 1 £(£+!) 

(Or*=nK--\nn + — M (2(1 + i)&/y/K) , st = „W= ( 15 7) 

4 2i y ' ' 2^-4/j, 

Notice that the first-order correction behaves as ln«, a fact which may be associated with gauge 
invariance. 

As with gravitational perturbations, the above analytic results are in good agreement with nu- 
merical results [29] (see ref. [27] for a detailed comparison). 

4 AdS/CFT correspondence and hydrodynamics 

A second unexpected connection comes from studies carried out using the Rela- 
tivistic Heavy Ion Collider, a particle accelerator at Brookhaven National Laboratory. 
This machine smashes together nuclei at high energy to produce a hot, strongly interact- 
ing plasma. Physicists have found that some of the properties of this plasma are better 
modeled (via duality) as a tiny black hole in a space with extra dimensions than as the 
expected clump of elementary particles in the usual four dimensions of spacetime. The 
prediction here is again not a sharp one, as the string model works much better than ex- 
pected. String-theory skeptics could take the point of view that it is just a mathematical 
spinoff. However, one of the repeated lessons of physics is unity - nature uses a small 
number of principles in diverse ways. And so the quantum gravity that is manifesting 
itself in dual form at Brookhaven is likely to be the same one that operates everywhere 
else in the universe. 

- Joe Polchinski 

There is a correspondence between JV = 4 Super Yang-Mills (SYM) theory in the large N 
limit and type-IIB string theory in AdSs x S 5 (AdS/CFT correspondence). In the low energy limit, 
string theory is reduced to classical supergravity and the AdS/CFT correspondence allows one to 
calculate all gauge field-theory correlation functions in the strong coupling limit leading to non- 
trivial predictions on the behavior of gauge theory fluids. For example, the entropy of ,jV = 4 SYM 
theory in the limit of large 't Hooft coupling is precisely 3/4 its value in the zero coupling limit. 

The long-distance (low-frequency) behavior of any interacting theory at finite temperature must 
be described by fluid mechanics (hydrodynamics). This leads to a universality in physical proper- 
ties because hydrodynamics implies very precise constraints on correlation functions of conserved 
currents and the stress-energy tensor. Their correlators are fixed once a few transport coefficients 
are known. 
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4.1 Hydrodynamics 

To study hydrodynamics of the gauge theory fluid, suppose it possesses a conserved current j^. For 
simplicity, let us set the chemical potential /J. = 0, so that in thermal equilibrium the charge density 
(j ) = 0- The retarded thermal Green function is given by 

G* v {<0,q) = -ijd 4 xe-^d(t)([j^x), j v (0)]>, (158) 

where q = (<a,q), x = (t,x). It determines the response of the system to a small external source 
coupled to the current. For small CO and \q\, the external perturbation varies slowly in space and 
time. Then a macroscopic hydrodynamic description for its evolution is possible [30]. 
For a charged density obeying the diffusion equation 

d Q f = DV 2 f, (159) 

where D is the diffusion constant with dimension of length, we obtain an overdamped mode with 
dispersion relation 

(0 = -iDq 2 , (160) 

The corresponding retarded Green function has a pole at co = —iDq 2 in the complex G)-plane. 

Another important conserved current is the stress-energy tensor T^ v . Its conservation law may 
be written as 



d T 00 + d i T 0i = o, 

d T 0i + djf ij = 0, 



(161) 



where 



f oo = r oo_ P; p = (T 0Q ) 
1 



P+P 



2 



pi = T ij - P 8 ij = tj d t T°J + djT 0i - 1 8 ij d k T Qk +t;S ij d k T 



3 



n()A' 



(162) 



and p (p) is the energy density (pressure) of the fluid, f] (£) is its shear (bulk) viscosity. 

One obtains two types of eigenmodes, the shear modes which consist of transverse fluctuations 
of the momentum density T , with a purely imaginary eigenvalue 

(0 = -iDq 2 , D = —^—, (163) 

P+P 

and a sound wave due to simultaneous fluctuations of the energy density T and the longitudinal 
component of momentum density T Ql , with dispersion relation 

CO = U s q-^-^- p ^ + 4 -^q 2 , U 2 = ^. (164) 

In a conformal field theory, the stress-energy tensor is traceless, so 

p = 3p , £=0 , u, = -j= (165) 
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4.2 Branes 

To understand the gravitational side of the AdS/CFT correspondence, consider a non-extremal 3- 
brane which is a solution of type-IIB low energy equations of motion. In the near-horizon limit 
r <^R where R is the AdS radius, the metric becomes 

ds 2 w = ± - (-f(u)dt 2 + dx 2 + dy 2 + dz 2 ) + 2 . , du 2 + R 2 dD. 2 5 , (166) 

li tyli J I li J 

2 

where T = -j^ is the Hawking temperature, and we have defined u = 4[, f{u) = \ — u 2 . The horizon 
corresponds to u = 1 whereas spatial infinity is at u = 0. 

According to the gauge theory/gravity correspondence, the above background metric with non- 
extremality parameter ro is dual to .yV = 4 SU(N) SYM at finite temperature T in the limit of 
N — > °o, gy M N — ► oo. For the retarded Green function 

G^ v ,x P ((0,q) = -iJd 4 xe-^ x d(t)([T^ v (x), T Ap (0)]). (167) 

we deduce by considering an appropriate perturbation of the background metric [30], 

Gxy^CO^q) = — {i2TlT(0 + q 2 ) . (168) 

leading to the shear viscosity of strongly coupled ,jV = 4 SYM plasma (Kubo formula) 

r\ = lim -L [dtdxe^'dTxyix), 7^(0)]) = ^N 2 T 3 . (169) 

»— >0 ZOO J a 

Other correlators may also be found by different perturbations of the metric. One obtains 



N 2 7lT 3 q 2 
%(i<D-®q* 

N 2 nT 3 co\q\ 
8(/co-^ 

N 2 7lT 3 CD 2 
SiiCO-^q 2 



G tXttx (CQ,q) - ,,.„ 7a&\ +---> 



G tx , xz {o),q) - - of:m ^-n, +- 



G xz ^{0),q) - 0/; „ ^, +---, (170a) 



where <& = -^. 



From the above results, one may deduce the viscosity tj. Indeed, recall from hydrodynamics 
Q> = —2-. Using the entropy 

3 n 2 t o 
s=-s Q = —N 2 T 3 , (171) 

where sq is the entropy at zero coupling, and the thermodynamic equations s = |^, p = 3p, we 
deduce p +p = ^-N 2 T 4 , therefore 

I, = ?-N 2 T 3 , 1 = i- (172) 

8 S 47T 

which agrees with the Kubo formula. It should be pointed out that there is no agreement unless 
s = js , a fact which is still poorly understood. 
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The above result for the viscosity is based on the gravity dual of the gauge theory fluid and 
should correspond to its strong coupling regime. At weak coupling, one obtains by a direct calcu- 
lation 

T7 1 

- » — (173) 

s 4n 

Thus the viscocity coefficient 77 varies as a function of the 't Hooft coupling, 

ri=fi 1 (g 2 Y M N)N 2 T 3 (174) 

where f^ (x) ~ _J lnx for x ^ 1 and /,, (x) = | for x S> 1 . 

4.3 Schwarzschild black holes 

In the metric considered above, the horizon was flat. This corresponds to the limit of a large black 
hole. For a black hole of finite size, the horizon is generally a sphere. Then the boundary of 
spacetime is S 3 x R. This may be conformally mapped onto a flat Minkowski space. Then by 
holographic renormalization, the AdSs-Schwarzschild black hole is dual to a spherical shell of 
plasma on the four-dimensional Minkowski space which first contracts and then expands (conformal 
soliton flow) [31]. 

Quasi-normal modes govern the properties of this plasma with long-lived modes (i.e., of small 
3 ft)) having the most influence. For example, one obtains the ratio 

v 2 1 (o 4 -40(o 2 + 12 . n(0 

-5- = —3i 5 — sin — - (175) 

8 6n co 3 -4co 2 



where V2 = (cos 20) evaluated at 6 = -| (mid-rapidity) and averaged with respect to the energy 
density at late times; 5 = : y 2 ~ J is the eccentricity at time t = 0. Numerically, ^ =0.37, which 
compares well with the result from RHIC data, % « 0.323 [32]. 

Another observable is the thermalization time which is found to be 

T = TRT7 ~ T^T- ~ °- 08 fm / c > r p-k = 300 MeV (176) 

not in agreement with the RHIC result % ~ 0.6 fm/c [33], but still encouragingly small. For com- 
parison, the corresponding result from perturbative QCD is % > 2.5 fm/c [34, 35]. 

The above results motivate the calculation of low-lying QNMs. Earlier, we calculated analyt- 
ically the asymptotic form of QNMs for large black holes. We obtained frequencies which were 
proportional to the horizon radius tq. We found an infinite spectrum, however we missed the lowest 
frequencies which are inversely proportional to r$. The latter are important in the understanding of 
the hydrodynamic behavior of the gauge theory fluid via the AdS/CFT correspondence. 

4.3.1 Vector perturbations 

We start with vector perturbations and work in a <f-dimensional Schwarzschild background. It is 
convenient to introduce the coordinate [36] 

u=(^Y~* (177) 



r 

The wave equation becomes 



-{d-?>) 2 u^f{u)(u^f{u)^'\ + V v {u)^> = cb 2x V , a = — (178) 

V / r 
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where prime denotes differentiation with respect to u and we have defined 



r- 



'0 

«, x I -o (d-2)(d-4) 2 . (« 

ri 4 



/W-^=l-.Af.-^) (179) 



W(m) = -£= f(u) { Z/ + ^ ^ '-u-**f{u) ^ °-^u } (180) 



where L 2 = -* — , — '- . 



First let us consider the large black hole limit tq — ► oo keeping ft) and L fixed (small). Factoring 
out the behavior at the horizon (u = 1) 

>F=(l- M )- i Af( M ) (181) 

the wave equation simplifies to 

s^F" + @ & F' + V A jF = (182) 

where 

■fi/ = -(d-3) 2 M« (1-M«) 



ld-% rf-1 , 

-(rf-3)[rf-4-(2rf-5)«3=5]ii3=5 -2(d-3) 2 



J— 1 1 — u 



n (d-2)[d-4-3(d-2)u^] 2 
L -\ u rf - 3 



^a,L - 4 

\2 a2 



2rf-8 , rf-1 , 



ft) , . _ s t ft) M <« (1 -M rf -3 



— + (^-3) 



1_ M ^ V ' (d- 1)2 (1- M )2 

jft) |d — 4— (2d — 5)w<*-3]w<<-3 2 i ffl u d-i (\ - ud-\ 



d-\ \-u v ' d-1 (l-«) 2 

We may solve this equation perturbatively by separating 

(J^, + ^)F = (183) 

where 

J%F = s4F" + S§ Q F' + %),oF 

JftF = (<% & -<% )F' + (V &L -tf 0fi )F 

Expanding the wavefunction perturbatively, 

F = Fo + F x + ... (184) 

at zeroth order the wave equation reads 

J%>F = (185) 

whose acceptable solution is 

F = u^ (186) 
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being regular at both the horizon (u = 1) and the boundary (u = 0, or *P ~ r 2 — ► as r — ► 00). 
The Wronskian is 

K««-3(l — W««-3 ) 



and another linearly independent solution is 



r 



F = Fol^ (188) 



^ 



which is unacceptable because it diverges at both the horizon (Fq ~ ln(l — u) for u rs 1) and the 



</-4 



boundary (F ~ m 2 ( d3 ) for u « 0, or *P ~ rT" — ► 00 as r — > 00). 
At first order the wave equation reads 

M'qFx = -JftF (189) 

whose solution may be written as 

fW fF J?{Fo 

The limits of the inner integral may be adjusted at will because this amounts to adding an arbitrary 
amount of the unacceptable solution. To ensure regularity at the horizon, choose one of the limits 
of integration at u = 1 rendering the integrand regular at the horizon. Then at the boundary (u = 0), 

F i=Fo - ' , - + regular terms (191) 

JO S&W 

The coefficient of the singularity ought to vanish, 

"^P=° < 192 > 

ssW 
which yields a constraint on the parameters (dispersion relation) 

a L 2 - mi d) - a 2 d) 2 = (193) 

After some algebra, we arrive at 

a = ^— ^ , ai=d-3 (194) 

(3 — 1 

The coefficient a 2 may also be found explicitly for each dimension d, but it cannot be written as a 
function of d in closed form. It does not contribute to the dispersion relation at lowest order. E.g., 
for d = 4, 5, we obtain, respectively 

65 1 5 1 

a2 = i08-3 ln3 ' 6"2 ln2 <195) 

Eq. (193) is quadratic in ft) and has two solutions, 

L 2 . d- 3 l? 

OX) « -i- — - , coi « -/ h i-j — 7 (196) 

a — 1 a2 a — 1 
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In terms of the frequency oo and the quantum number £, 

.£{£ + d-3) o)i 



<M[) 



(d-l)r 



.d-3 .£(£ + d-3) 



n) 



a 2 



(d- 



iH 



(197) 



The smaller of the two, ooo, is inversely proportional to the radius of the horizon and is not included 
in the asymptotic spectrum. The other solution, C0\ , is a crude estimate of the first overtone in the 
asymptotic spectrum, nevertheless it shares two important features with the asymptotic spectrum: 
it is proportional to tq and its dependence on £ is G{\jr^). The approximation may be improved 
by including higher-order terms. This increases the degree of the polynomial in the dispersion 
relation (193) whose roots then yield approximate values of more QNMs. This method reproduces 
the asymptotic spectrum derived earlier albeit not in an efficient way. 

To include finite size effects, we shall use perturbation theory (assuming 1/ro is small) and 
replace <ffi[ by 

J^J = M\ H — jJ^; (198) 



where 



,^ H F = £/ H F" + S8 H F' + VhF 



(199) 



The coefficients may be easily deduced by collecting &{\jr^) terms in the exact wave equation. 
We obtain 






-2{d-3) 2 u 2 {\-u 
-(d — 3)u 



(d - 3) (2 - 3m) -(d- 1) L-w^ 



d-l 
1 —U<l-3 



% 



d-l 



H 



d-A-{2d-5)u- {d-l) — W" 3 

1-m~ 



Interestingly, the zeroth order wavef unction Fq is an eigenf unction of J^k, 

jehFo = -(d-2)F 
therefore the first-order finite-size effect is a simple shift of the angular momentum operator 

s^ 'y s^ o CI _ 



,2 
r () 



The QNMs of lowest frequency are modified to 

.£{£ + d-3)-{d-2) 



(Oo 



For d = 4, 5, we have respectively, 



(Oo 



(d-l)ro 



3r 



+ ^(lAo) 



_.('+!) 



4r 



(200) 



(201) 



(202) 



(203) 



in agreement with numerical results [29, 31]. 

According to the AdS/CFT correspondence, dual to the AdS Schwarzschild black hole is a 



gauge theory fluid on the boundary of AdS (5' 



■d-l 



). Consider the fluid dynamics ansatz 



-iflTv 



J^e' lL "Yi 



(204) 
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where w,- is the (small) velocity of a point in the fluid, and V,- a vector harmonic on S d ~ 2 . De- 
manding that this ansatz satisfy the standard equations of linearized hydrodynamics, one arrives at 
a constraint on the frequency of the perturbation Q. which yields [37] 

»--' <(<+ 1;V'- 2 W '3> P05) 

{a— l)r 
in perfect agreement with its dual counterpart. 

4.3.2 Scalar perturbations 

Next we consider scalar perturbations which are calculationally more involved but phenomenologi- 
cally more important because theirspectrum contains the lowest frequencies. For a scalar perturba- 
tion we ought to replace the potential Vy by 

-2 



Vs(«) 



/(«) 



m+ ( 1 + — ) u 



1 

7 2 

{. / 1 \ 2d-8 , , , / 1 \ d-5 

d(d - 2) ( 1 + -j J u^ - 6(d -2)(d - 4)m l\ + — ] u~^ 

+ (d-4)(d-6)m 2 u-^ + (d-2) 2 (\ + \j u 3 



1\ 2 

r 2 

'o 



+2(2d 2 -lW+18)m(l + — ) z/ 2 

(d-4)(d-6)(\ + j,) 2 / n 

+ j-^ ^ M 2 -3(d-2)(d-6)m 2 1 + 3 U 

6(</-2)(</-4)m(l + £) ^2 

. i °-^u + 2(d-l)(d-2)m 3 +d(d-2)^ r 



where m - o^l^ti^ 2 ! - 2(l+<*-2)(l-i) 
where m " 2 (d-l)(rf-2)^ " {d-X){d-2)4 ■ 

In the large black hole limit ro — ► °° with m fixed (small), the potential simplifies to 



(206) 



d-\ 



?A°)i 



1 — «d-3 



''_5 
-3 



b 4(m + u) 2 1 

+( J - 4) ( J - 6)m V ^ + {d- 2) V 

+2(2d 2 - 1 Id + 18)/ro 2 - 3(d - 2)(d - 6)m 2 u + 2(d -!)(</- 2)m 3 



(207) 



The wave equation has an additional singularity due to the double pole of the scalar potential at 
u = —m. It is desirable to factor out the behavior not only at the horizon, but also at the boundary 
and the pole of the scalar potential, 



^I-^t-FW (208) 

m + u 

28 



Then the wave equation reads 
where 

. o 2rf-8 , d-l „ 

^ = _(j_3) 2 M in-(i_ M 3^; 

m & = -(d-3)u-3=*(i-uz=i) 



£fF" + & & F' + <tf & F = 



(209) 



d-4 2(d-3) 



m + u 



d-l, d-5 
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4u 2 



u(m + u) [m + u) 
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2d-& d-l , 
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(J-l) 2 (1-w) 2 
We shall define zeroth-order wave equation as J%Fq = 0, where 

The acceptable zeroth-order solution is 

Fa(u) = 1 



(210) 



(211) 



which is plainly regular at all singular points (u = 0,1,— in). It corresponds to a wavefunction 

d-4 

vanishing at the boundary QV ~ r 2 as r — ► 00). 
The Wronskian is 



>T: 



(m + w) 



2rf-8 rf-1 , 

W<*-3 (1 -W-3 



(212) 



and an nacceptable solution is Fo = f #^. It can be written in terms of hypergeometric functions. 
For d > 6, it has a singularity at the boundary, Fq ~ u~^ for u « 0, or 1* ~ r~ — ► 00 as r — > 00. 
For d = 5, the acceptable wavefunction behaves as r -1 / 2 whereas the unacceptable one behaves as 
r _1 ' 2 lnr. For d = 4, the roles of Fq and Fq are reversed, however the results still valid because the 
correct boundary condition at the boundary is a Robin boundary condition [36, 37]. Finally, we note 
that Fq is also singular (logarithmically) at the horizon (u = 1). 

Working as in the case of vector modes, we arrive at the first-order constraint 

- 1 % 



because JfiFi 



1^0 



/o £?W ° 

^0)^0 + %-^b = %• This leads to the dispersion relation 



ao — aiJG) — a2&) 
After some algebra, we obtain 

d-l l + (d-2)m 



S.2 







ao 



(1+m 



$A2 



ai 



<i — 3 1 , 

, a 2 = -{l + 0(m)} 



(1+m) 



m 



(213) 



(214) 



(215) 



29 



For small m, the quadratic equation has solutions 



„ + d — 3 /d— 1 „ 

^o ~ -«— 2~ m ±y — 2~ m ( 216 ) 

related to each other by 6)q = —0)q*, which is a general symmetry of the spectrum. 

Finite size effects at first order amount to a shift of the coefficient ao in the dispersion relation 

a ^a + ^a// (217) 

After some tedious but straightforward algebra, we obtain 

a H = -{l + 0{m)} (218) 

m 



The modified dispersion relation yields the modes 



„_i_ d — 3 Id — 1 

&Q*t-i-;-m±J-z-m + \ (219) 



In terms of the quantum number (., 



... ^£(l + d-3)-(d-2) £(£ + d-3) 



<«-m-3) V-w-2)h. ^^^ <220) 

in agreement with numerical results [31]. 

Notice that the imaginary part is inversely proportional to ro, as in vector case. In the scalar 
case, we also obtained a finite real part independent of ro. It yields the speed of sound v s = , l 
which is the correct value in the presence of conformal invariance. 

Turning to the implications of the above results for the AdS/CFT correspondence, we may 
perturb the gauge theory fluid on the boundary of AdS (S d ~ 2 x W) using the ansatz 

Ui = J^e- iClr Vi$ , Sp = J^'e- ia *$ (221) 

where m is the (small) velocity of a point in the fluid and 8p is a pressure perturbation. They are 
both given in terms of S, a scalar harmonic on S d ~ 2 . Demanding that this ansatz satisfy the equations 
of linearized hydrodynamics, one obtains a frequency of perturbation Q. in perfect agreement with 
our analytic result [36, 37]. 

4.3.3 Tensor perturbations 

Finally, for completeness we discuss the case of tensor perturbations. Unlike the other two cases, 
the asymptotic spectrum of tensor perturbations is the entire spectrum. To see this, note that in the 
large black hole limit, the wave equation reads 

-{d-3) 2 {u^ -u 3 )^" -{d-3)[{d-A)u^-{2d-5)u 2 ]^>' 

n d{d-2) _ A (d-2) 2 cb 2 > n) 
+ <L + if ^-3 + if jzt^ 

* * 1— W-3 



30 



For the zeroth-order equation, we may set L = = ft). The resulting equation may be solved exactly. 
Two linaerly independent solutions are QV = Fq at zeroth order) 



F (u)=uW=V , F (u) = u 53=3) in M _ M a=3 ) (222) 



d-\ 



Neither behaves nicely at both ends (u = 0, 1). Therefore both are unacceptable which makes it 
impossible to build a perturbation theory to calculate small frequencies which are inversely propor- 
tional to ro. This negative result is in agreement with numerical results [29, 31] and in accordance 
with the AdS/CFT correspondence. Indeed, there is no ansatz that can be built from tensor spherical 
harmonics T,y satisfying the linearized hydrodynamic equations, because of the conservation and 
tracelessness properties of T,y. 

5 Conclusion 

We discussed the calculation of analytic asymptotic expressions for quasi-normal modes of various 
perturbations of black holes in asymptotically flat as well as anti-de Sitter spaces. We also showed 
how perturbative corrections to the asymptotic expressions can be systematically calculated. 

In view of the AdS/CFT correspondence, in AdS spaces we concentrated on low frequency 
modes because they govern the hydrodynamic behavior of the gauge theory fluid which is dual to 
the black hole. Thus, these modes provide a powerful tool in understanding the hydrodynamics of 
a gauge theory at strong coupling. They may lead to experimental consequences pertaining to the 
quark-gluon plasma produced in heavy ion collisions at RHIC and the LHC. 
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